Coulomb Interactions and Nanoscale Electronic Inhomogeneities in Manganites 
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We address the issue of endemic electronic inhomogeneities in manganites using extensive sim- 
ulations on a new model with Coulomb interactions amongst two electronic fluids, one localized 
(polaronic), the other extended (band-like), and dopant ions. The long range Coulomb interactions 
frustrate phase separation induced by the strong on site repulsion between the fluids. A single 
quantum phase ensues which is intrinsically and strongly inhomogeneous at a nano-scale, but ho- 
mogeneous on meso-scales, with many characteristics (including colossal responses)that agree with 
experiments. This, we argue, is the origin of nanoscale inhomogeneities in manganites, rather than 
phase competition and disorder related effects as often proposed. 

PACS numbers: 71.10.-w, 71.27.-|-a, 75.47.Lx 



An intriguing generic characteristic of complex solids 
such as doped manganites 0, 0, I3i cupratesQ and 
cobaltates|5| seems to be the coexistence of patches 
of metallic and insulating regions (often dubbed as 
'phases'). This 'electronic inhomogeneity' can vary from 
nanometers to microns and can be static or dynamic P, 
IE IE ^ 111' Questions as to whether this can be character- 
ized as 'electronic softness' and is a defining feature of 
these materials, its origins and role in determining their 
electronic properties, etc., are among the most actively 
explored issues. However, proposed mechanisms |3,[3 ^^''^ 
either at the level of scenarios or toy models, and do not 
adequately address the specifics of the scale and nature 
of inhomogeneities in the actual systems. 

Motivated by the fact that electronic inhomogeneities 
came to prominent attention first in manganites 0, in 
this letter we address these issues using extensive simu- 
lations on a new, "extended £b model". This is an ex- 
tension of the £b model, involving two kinds of electronic 
fluids, one localized and polaronic {£), the other extended 
and band-like (6), which was introduced recently and 
shown to describe successfully many puzzling phenomena 
(including colossal magnetoresistance) observed in doped 
manganites of the type Rei-a^Ak^^MnOs (Re=rare earth 
ions such as La, Nd; Ak= alkaline earth ions, such as 
Ca, Sr). The extensions invoke long range Coulomb in- 
teractions amongst the two fluids, and dopant ions. We 
explore here their consequences, especially in the con- 
text of issues connected with electronic inhomogeneities 
where they play a crucial role. 

The ib model|E,ll] is based on the idea that under the 
conditions prevailing in the doped manganites, the elec- 
trons populating the doubly degenerate eg states cen- 
tered at the Mn sites spontaneously reorganize them- 
selves into two types of electron fluids which coexist. One 
is obtained by populating essentially site localized states 
labeled £ which are polaronic, with strong local Jahn- 
Teller(JT) distortions of the oxygen octahedra, an energy 
gain EjT ~ 0.5 eV[lOj and exponentially reduced inter- 



site hopping. The other, labeled b, is a fluid of broad 
band, non-polaronic electrons, with no associated lattice 
distortions, and undiminished hopping amplitudes. For 
a generic x in the regime (0.1 < x < 0.4) of interest 
to us in this paper, manganites do not exhibit orbital 
long range order and can be regarded as 'orbital liquids'. 
Hence one can characterize the hopping of the b elec- 
trons by a single orbitally averaged number t 0.2) 
eV and ignore the eg orbital index. There is a strong 
local repulsion between the two fluids, as double occu- 
pancy on a polaronic site costs an extra Coulomb energy 
U (~ 5eV) . The spins of £ and b are enslaved to the Mn- 
t2g spins (5 = 3/2) on site due to the large ferromagnetic 
Hund's coupling Jh (~ 2eV). Furthermore, there is a 
new, occupancy dependmt, ferromagnetic nearest neigh- 
bor exchange coupling^ between the t2g core spins, of 
order x{l — x)t^ /Ejt, referred to as "virtual double ex- 
change" ; which overwhelms the normal super-exchange 
for X > 0.1. Hence, as is seen in manganites over the 
above mentioned range of doping x, the ground state is 
ferromagnetic (insulating or metallic). In the simplest 
picture, assuming all the t2g spins and the Cg spins to be 
aligned parallel, the £ and b electrons can be regarded 
as spin-less, leading naturally to the Falicov-Kimball 1 1 ij 
like £b Hamiltonian 
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Here f f and 6] create £ polarons and b electrons respec- 
tively at the sites i of a cubic Mn lattice, and na = £\£i 
and nu = b\bi are the corresponding number operators. 

In this Letter we extend the above model to include 
the long range Coulomb interactions that are necessarily 
present in the doped manganites. This is done at the 
simplest level by associating quenched (spatially fixed) 
charges — |e| (with respect to the Re^"*" background) at 
a fraction x of random, body centered, 'Ak sites', and 
annealed charges -|-|e| with the corresponding deficits in 
the Mn-Bg electron occupancy (overall reduced to (1 — x) 

per site). In terms of the hole operator [h\ = £i which 
removes an £ polaron at site i), and the electron charge 



2 




FIG. 1: (color online) Real space electronic distribution ob- 
tained from simulations on a 16^ cube. Magenta (darkest) de- 
notes hole clumps with occupied b electrons, white (lightest) 
denotes hole clumps with no b electrons, cyan (2nd lightest) 
denote singleton holes, and light blue (2nd darkest) repre- 
sents regions with i polarons. The configuration on the left 
shows isolated clumps with occupied 6-electrons (ti-electron 
puddles). For larger doping, percolating clumps are obtained 
and the system is a metal (right). The inset is results in the 
absence of long range Coulomb interaction (Vo = 0.0) and 
shows 'macroscopic phase separation'. All energy scales are 
in the units of t. 




FIG. 2: (color online) Critical doping levels Xd and Xc2- Xci 
separates the lightest region, where no b states are occupied, 
from the intermediate shade region, which has b states occu- 
pied in puddles. The darkest shaded region, with x > Xc2, has 
some occupied b states that percolate through the 10^ box. 
The solid line corresponds to the DMFT result[^ for Xci- 



operator qi = h\hi — b\bi, which has the average value x 
per site because of overall charge neutrality, the extended 
model Hamiltonian is 

H = m, + Hc,Hc^T^.q. + ^T^. (2) 

The Coulomb term He has two parts; the charge at site 
i has energy qi^i, where $i is the electrostatic potential 
there due to Ak^+ ions, and the interaction between the 
charges at site i and j leads to an energy Vq . 

We now describe the results of numerical simulations of 
the above model Q on finite 3d periodic lattices of sizes 
up 20'^. The electrostatic interaction is treated in the 
Hartree approximation, using the mean b electron charge 
at site i , i. e., qi is approximated by (qi) = h\hi — {b\bi). 



Furthermore, since U is large, for simplicity we take the 
U oo limit; i.e., assume that b electrons do not hop 
to sites occupied by £, so the two kinds of electrons form 
disjoint clusters. On a hole-cluster, which has two or 
more /i-sites each accessible to the other by electron hops, 
which we will refer to as a "clump", the h electron en- 
ergy levels are determined exactly by solving the intersite 
Hamiltonian iff, = —tY^ij{b\bj+h.. c). The ground state 
(i.e, at zero temperature) is obtained by starting from an 
overall charge neutral trial configuration random of £ po- 
laronic sites, and performing electron transfers that lower 
energy, till none such exist. Some of these are of the type 
£ — /i as in a classical Coulomb glass ^3 j which an ^ 
electron is moved to a hole site. In addition there are 
£ — b and b — b transfers involving one or two 6-states 
whose energy is quantum mechanically obtained for a 
particular clump structure. At each iteration then, we 
find the best possible transfer, i. e., the one which low- 
ers most the (occupied) single particle levels whose ener- 
gies include Coulomb interactions. Then we perform the 
transfer, and update the clump structure, if necessary, 
i. e., redistribute the b electrons in the new clumps. The 
process is repeated till the final 5-clump/£-polaron struc- 
ture is stable against all further electron transfers. This 
is a new generalization of the common Coulomb glass 
simulation]!^] which includes the quantum mechanically 
obtained b states within their clump or puddle. The elec- 
trostatic energy is calculated accurately using the Ewald 
technique and fast Fourier transform routines. 

In the discussion below, all length scales are normal- 
ized by the lattice parameter a, and we use dimension- 
less energy parameters Ejt and scaled by the hop- 
ping amplitude t. Realistic values jlOf for manganites are 
2.0 < Ejt < 3.0 and 0.01 < Vq < 0.1 with t w 0.2eV. 
In the absence of Vq, the system "phase separates", i. e., 
holes move to one side of the box and several £ polarons 
are converted to b electrons that occupy the low energy 
band states with energies below —Ejt in this large clump 
(FIG. ni see inset). This phase separation is due to 
strong local correlations (large U between £ s and b s) 
and is in agreement with known results for the Falicov- 
Kimball model[lll|. The presence of Vq renders this phase 
separation energetically unfavorable, and intermixes the 
phases as expected The favored electronic configu- 
rations strongly depend on the JT energy Ejt and the 
doping X. Two examples are shown in FIG. ^ For a 
given Ejt, there are two critical values of doping, Xd 
and Xc2, as shown in FIG. 12 For x < Xci, there are 
no occupied b-electron states in the system, and the holes 
form a Coulomb glass^. For doping larger than Xci, 
some occupied 6-electron puddles appear (FIG.Q^a)). On 
further increase of doping, clumps interconnect and per- 
colate through the simulation box (FIG. ^ (b)), giving 
rise to the possibility of occupancy of 6-states extended 
throughout the system, and hence metallicity. The re- 
sults for Xci and Xc2 are insensitive to Vq for Vq < 0.5 
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FIG. 3: (color online) Density of states (DOS) of £ polarons 
and b electrons (average over 100 random initial conditions, 
size: 10"^). The occupied states are shaded. The chemical 
potential fi and the b band center Ei, are marked. Inset: soft 
Coulomb gap of the polarons. 

(the typical value of 0.01 < Vq < 0.1 in manganites falls 
in this range). For larger values of Vq, the clumps are 
always too small, and are never occupied by electrons, 
and the system is an insulator for all x < 0.5 . 

The physics behind these results is uncovered by a 
study of the density of statcs(DOS) of the ^-polarons and 
6-electrons shown in FIG. |2| Wc find that the chemical 
potential fi (FIG. 13 goes essentially as ^ = —Ejt + Vq. 
The polarons form a Coulomb glass which has a soft 
gap (FIG. 13 inset) at the chemical potential, as in the 
classical Coulomb glass [l^. Furthermore, the polarons 
have their energies spread over an energy scale of or- 
der Vb.[i3 The b electrons have a band-like DOS with a 
bandwidth that depends on the doping x. In fact, we find 
that the effective half bandwidth of the h electrons scales 
as Deff = -Dov^ (^0, the bare half-bandwidth, = Qt 
for the cubic lattice), in remarkable agreement with an 
earlierj^ single site dynamical mean field thcory(DMFT) 
result obtained for (^! Furthermore, the simulation re- 
sults show that the h band center (marked by in the 
figure) is not affected by doping, and scales as Ef, « Vq 
(for small Vb < 0.5). 

The behavior of the chemical potential and the band 
center are due to the fact that the random distribution of 
Ak ions plays an important role in determining the real 
space structure. From a study of the positional correla- 
tion function between two holes, and a hole and an Ak 
ion, we find that the holes tend to cluster around the Ak 
ions^3- Furthermore, the electrostatic screening is quite 
strong in this system in that these correlation functions 
(both hole-hole and hole-Ak ion) reach a plateau within 
a few lattice spacings. Since the holes cluster around the 
(oppositely charged) Ak ions, it is clear that a typical 
(. polaron will have a larger number of £-polaron neigh- 
bors than hole sites. Thus, the t polaron sites see, on an 
average, an electrostatic potential of order Vq, and the 
average energy of an i polaron is increased from —Ejt 



to —Ejt + Vq. Since the chemical potential of a Coulomb 
glass is the average energy of the states ^3 1 it follows that 
the chemical potential is —Ejt + Vq. Similarly every hole 
site sees an electrostatic potential of Vq, implying that 
the b band center is placed at Vq. These observations 
and the fact that the b bandwidth scales as ^/x, suggests 
that the minimum doping for the the 5-band bottom to 
touch the chemical potential (i.e., for 6-state occupancy) 

is Xci = ) , exactly the DMFT prediction[|| for ||TJ! 

However, we note that the simulation results include 
important new physics not contained in the DMFt|^, 
namely that b state occupancy does not by itself make 
the system a metal, as these states are localized inside 
the clumps. Based on inverse participation ratio, geo- 
metric percolation of the clumps and Kubo conductivity 
calculations [il^. we have good estimates for a second, 
higher, critical doping level, Xc2, at which the system ac- 
tually becomes a metal (FIG.EJ, for which the occupied 
b-states should extend across the simulation box. Our re- 
sults suggest that for x < Xci , electrical transport should 
be that of a classical Coulomb glass (tr ~ e~^^^). For 
Xci < x < Xc2, the transport should have two contri- 
butions - the first a Coulomb glass contribution of the 
£ polarons and the second a variable range type inter- 
puddle hopping of the electrons in the b puddles, as is 
indeed observed in doped manganites|l3j . For x > Xc2 
we find a highly resistive metal. 

Another interesting aspect that we have 
investigated j 1 8^ is the clump size R, and its depen- 
dence on Vb- For a (fictitious) uniform distribution of 
Ak ions, one can show by an approximate analytical 
calculation that R ^ "^/^' corresponding, for Ejt = 2.5 
and X = 0.3, to clump sizes between 10 and 5 lattice 
spacings for Vq between 0.01 and 0.1. The clump sizes 
for the more realistic, random, distribution of Ak ions 
obtained from our simulations are much smaller; even 
a very small Vq produces clumps that are four to five 
lattice spacings, and this spacing is essentially indepen- 
dent of Vq for realistic values of Vq - in stark contrast 
to the analytical result above. Thus the long range 
Coulomb interaction is a 'singular perturbation' that 
prevents macroscopic phase separation; but the sizes 
and the distribution of the clumps are determined by the 
random distribution of the Ak ions. Thus we conclude 
that doped manganites as modeled by Q are necessarily 
and intrinsically electronically inhomogencous, on a 
nanometric scale. 

Our results have been obtained for a Hamiltonian and 
energy parameters that are very realistic especially in 
low bandwidth manganites with a large ferromagnetic 
region in their phase diagram. In sharp contrast to 
some of the earlier scenarios proposed in manganites, 
the nanoscale inhomogeneities we obtain are not due 
to 'phase competition' induced 'phase separation' be- 
tween 'insulating' and 'metallic' phases frustrated by dis- 
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order as suggested by studies on 'toy models', such as 
spin Hamiltonians0 or Hamiltonians with two locaHzed 
states and electron lattice coupHngQ- Rather, they arise 
due to the the long ranged Coulomb interactions frustrat- 
ing the phase separation induced by strong local correla- 
tions. The mechanism itself has been discussed in a vari- 
ety of contexts 0, Q, but to the best of our knowledge, 
ours is the first quantitative study on a realistic model 
(that includes dopant ions) for any correlated oxide. 

We emphasize that the nanoscale inhomogeneities we 
obtain are present in both the insulating and metallic 
phases of doped manganites between which one has a 
transition as a function of doping, at Xc2', and further- 
more, as is clear from Figs. H and El each of these con- 
stitutes a valid, single, thermodynamic phase that is ho- 
mogeneous on meso-scales. These results are in confor- 
mity with the electron holography results of Loudon et 
al. [l5j , where even the ferro-metallic state is seen to have 
interspersed in it nanoscale 'insulating regions', which 
in the context of the present simulation are just the £- 
clusters |2l|. 

Furthermore, our work suggests that mesoscale phase 
separation^], or proximity to multicriticalitv \l^ ] . are not 
essential for explaining CMR in manganites. Given the 
correspondence between the simulation results for the 
extended ih model and the (homogeneous) DMFT re- 
sults for the simple model, the single metallic phase 
obtained here will show a ferro-metal to para-insulator 
transition as well as CMR due to strong Hund's cou- 
pling between the Mn <2g core spins and the b elec- 
trons and thermal fluctuations. As shown elsewhere 0, 
the CMR arises from small field induced changes in the 
energetics of the extended b states which are responsi- 
ble for charge transport, and the consequent exponential 
changes in their occupancy. Indeed, by way of experimen- 
tal confirmation, we note that there are many mangan- 
ites without mesoscale inhomogeneities that show colos- 
sal responses 1 16j. Finally, our work reinforces sugges- 
tions that the mesoscale patterns seen in experiments 
arise from other sources such as long range elastic strains, 
possibly due to defects^, 0|. In our model, they can be 
generated from the strain dependence of the local energy 
parameters such as Ejt and t. Investigations along these 
lines are in progress. 
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